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Abstract 



This work closes certain gaps in the hterature on material reference systems in 
general relativity. It is shown that perfect fluids are a special case of DeWitt's 
relativistic elastic media and that the velocity-potential formalism for perfect 
fluids can be interpreted as describing a perfect fluid coupled to a fleet of 
clocks. A Hamiltonian analysis of the elastic media with clocks is carried 
out and the constraints that arise when the system is coupled to gravity are 
studied. When the Hamiltonian constraint is resolved with respect to the clock 
momentum, the resulting true Hamiltonian is found to be a functional only of 
the gravitational variables. The true Hamiltonian is explicitly displayed when 
the medium is dust, and is shown to depend on the detailed construction of 
the clocks. 
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I. INTRODUCTION 



The use of material reference systems in general relativity has a long and noble history. 
Beginning with the systems of rods and clocks conceived by Einstein and Hilbert 0, 
material systems have been used as a physical means of specifying events in spacetime and 
for addressing conceptual questions in classical gravity. That such systems also provide 
important tools for quantum gravity was pointed out by DeWitt [Q, who used them to 
analyze the implications of the uncertainty principle for measurements of the gravitational 
field. 

The original systems of rigid rods and massless clocks discussed by Einstein and Hilbert 
represent unphysical idealizations. Since their time, attempts have been made to remedy this 
shortcoming by developing a more physically realistic description of the reference medium. 
While still a phenomenological description, a dynamical reference system consistent with 
relativity may be found in the clocks and elastic media studied by DeWitt ||^. In the same 
spirit, perfect fluids have been employed as reference systems in the quantization of 
various model problems in gravity |^. More recently, Kuchar et al. P] have developed a 
scheme for incorporating reference systems in general relativity through the introduction 
of coordinate conditions For the cases examined thus far, the reference materials that 
arise through this approach have certain unphysical properties. This has motivated the 
investigation of pressureless perfect fluid (dust), which is unrelated to any obvious coordinate 
condition, as a phenomenological but physically realistic reference system [|]. 

Here we pursue the modest goal of closing certain gaps in the literature on reference 
materials. First, we establish in Sec. II the connection between the elastic media of DeWitt 
and the relativistic perfect fluids of Refs. 0,^,0], which include dust as a special case. 
In particular we show that perfect fluids are equivalent to elastic media when the later are 
homogeneous and isotropic. 

We then turn to the study of reference systems defined by the coupling of clocks to 
reference materials. In Sec. HI. A we demonstrate that a perfect fluid, as described by the 
action S of Refs. ||TO|JTl|l for the isentropic case, can be interpreted as a perfect fluid coupled to 
a fleet of clocks (with the details presented in Appendix A). We show the inequivalence of this 
clock coupling to the coupling used by DeWitt |^, and discuss the advantages of DeWitt 's 
method. We also point out that for a nonisentropic perfect fluid S, the thermasy |T^ can 
be reinterpreted as a clock variable with DeWitt-type coupling. Then, in preparation for 
Sec. IV, we perform in Sec. III.B a Hamiltonian analysis of the reference system formed by 
adding clocks (in the manner of DeWitt) to the elastic medium. This proceeds along the 
lines of that done in Refs. [T^JllI for perfect fluids. 

In Sec. IV we study the coupling of the reference system to gravity and the resulting 
canonical constraints. When the Hamiltonian constraint is resolved with respect to the clock 
momenta, we find that the true Hamiltonian depends only on the gravitational variables, 
not on the clock or particle variables. (Some details are derived in Appendix B.) The case 
of dust is particularly simple, and we explore the use of more general clocks than those 
used in Ref. [||. For this medium and for appropriately designed clocks, the constraints 
can be resolved with respect to the clock momenta using only analytic operations (i.e., 
without taking a square root). This allows us to bypass a technical difficulty when defining 
the quantum theory § although, as we point out, interpretational difficulties appear in its 
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place. 

We use the following notation, which is consistent with that of Refs. [10,11 1. The action 
for a reference system, which includes both matter and clocks, will be denoted by 5", while 
the action for the corresponding reference material alone, without clocks, will be denoted 
by S. Thus, the addition of clocks changes a 'barred' system into an 'un-barred' system. 



II. ELASTIC MEDIA AND PERFECT FLUIDS 

In this section we address the relation between De Witt's elastic media and the perfect 
fluids of Refs. p|-[ll|. In fact, the main difference is that the elastic medium is presented 



in the Lagrangian picture while the perfect fluids are typically presented in the Eulerian 
picture. Thus, we first review the elastic medium in the original Lagrangian description, 
and then rewrite it in the Eulerian picture. The perfect fluid action is then recognized to 
be a special case of the action for an elastic medium. 

We use the terms Lagrangian and Eulerian in the following way. The term "Lagrangian 
picture" refers to the description in which the basic variables tell the spatial location of 
a given particle, or the spacetime location of a given event on a given particle world line. 
The term "Eulerian picture" refers to the description in which the basic variables tell which 
particle resides at a given spatial location, or which particle passes through a given spacetime 
event. 

A. A Review of Elastic Media 

We now begin with a summary of DeWitt's elastic media 0. A single free relativistic 
particle moving in a spacetime Ai = H x M with metric 7q,/3 can be described by the action 

^i[T";7a/3] = - / t/ormv/-t"t'37„^(T) . (2.1) 

Here, a is an arbitrary parameter along the particle world line and y°' = T°((t) are the 
spacetime coordinates of the particle. Also, m > is the mass of the particle and the dot 
denotes a derivative with respect to a. The semicolon notation in ^ifT"; 7q,^] indicates that 
this action is to be varied with respect to T", with •ja/s treated as a background field. 

As in Ref. 0, we may also consider fleets of such particles and we may add local inter- 
actions. If the particles are labeled by a set of Lagrangian coordinates (\ i G {1, 2, 3}, then 
such a system can be described by the action 

^..[T";7a/3] = / da jj'c[Hmi + w)4^¥¥^^^)^ , (2.2) 



where S is the "matter space" manifold [§|10| whose points C, & S label the particle world- 
lines. In Eq. (2.2), the dynamical variables T" are functions of a and The quantity n 
is the particle number density, so that nd^C, is the number of particles in the coordinate 
cell d^C,. The quantity w_ is the interaction energy density in the co-moving frame, so that 
w_d^( is the interaction energy in the coordinate cell d^( as measured in the rest frame of 
the particles. In order to make the transformation properties of these expressions clear, we 
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explicitly indicate density weights with respect to changes of the coordinates C ^ with 
underlines. Thus, n and w are densities in the matter space S. 

The functions m, n, and w can depend explicitly on and, in addition, the interaction 
energy w can have an ultralocal dependence on the matter space metric ("fleet metric") hij. 

The fleet metric is defined such that ds = ^ hijdQd(^^ measures the orthogonal distance ds 
between neighboring world lines with Lagrangian coordinates C C + dC- By including 
an explicit dependence on Q in the mass m and interaction energy w, we allow for the 
possibility that the particles are not identical. In terms of the matter four-velocity 

U"" = I . . , (2.3) 

the matter space metric takes the form 

h^, = T^, (7„^ + f/«f//3)T^,- (2.4) 

where the commas denote derivatives with respect to In general hij will not be the metric 
of any hypersurface of the spacetime manifold M.. 

The system described by the action S'^^ of Eq. ( |2.2| ) is referred to as an elastic medium. 
The subscript LL indicates that it is the Lagrangian form (as opposed to the Hamiltonian 
form) of the action in the Lagrangian picture. Observe that this system is reparametrization 
invariant. That is, S^l is invariant under the transformation = — T"e induced by 
a reparametrization a — >■ a + e(cr, Q of the particle worldlines, where e vanishes at the 
endpoints in cr. 



B. The Connection to Fluids 

We will now derive the Eulerian description of the elastic medium, in which the action is 
written as an integral over arbitrary spacetime coordinates on M., and show its relations 
to the perfect fluids of Refs. |P-[TT||. 

To begin, let us assume that in the region of spacetime described, one and only one 
particle of the medium passes through each event. Then the Lagrangian coordinates Q 
along with the worldline parameters a form a set of coordinates in the spacetime region. 
These coordinates are related to the coordinates by the mappings = T"(cr, Q. We can 
introduce the inverse mappings 

<y = Z\y), C = Z\y), (2.5) 

such that = T°(Z°(?/), Z*(?/)) is an identity. Note that Z'^{y) give the labels of the 
particle present at the event y". The transition from the Lagrangian picture to the Eulerian 
picture is obtained by a change of dynamical variables in which T"(ct, Q is replaced by Z^{y) 
and Z\y). 

To perform this change of variables, we first calculate the Jacobian \d{a,()/dy\. Let 
us assume that the coordinate system a, has the orientation of Ai. Then we have the 
identity 
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da A dC^ A dC^ A dC^ = Z\f3 Z^,^ Z^s dy"^ A rf/ A dy'^ A c/y^ 



1 ^.,-PiS^^^^2'^^ ZJ^^ Z\, dy"^ A dy^ A rfy^ A dy"" , (2.6) 



3! 

where the commas followed by greek letters denote derivatives with respect to y". Here, 
^a/B-yS jg ^YiQ totally antisymmetric contravariant tensor on A4 with e^-"^^^ = —1/^—7, and 7 
is the determinant of the spacetime metric jais- Similarly, eijk is the totally antisymmetric 
covariant tensor on S with £123 = Vh. Observe that the inverse fleet metric can be written 
as 

h'^ = Z\^j^fZ^,p , (2.7) 

so that h can be expressed in terms of the Eulerian variables as h = 1/ det(/i*-' ). It will also 
be convenient to express the particle four-velocity f/" in terms of the new variables: 

U'' = -(l/3!)e"^^^Z\^ Z\., Z\s e^,k ■ (2.8) 

One can see that this expression is indeed the particle velocity by verifying that U^Z^,^ = 
and V^Un = — 1- This allows us to rewrite the measure (12.61) as 



da A dC} A dC"^ A dC,^ = ^j-^/h f/°rf|/° A dy^ A dy^ A dy^ . (2.9) 

Finally, note that t" is proportional to f/" and t"Z°,„ = 1, so that t" = f/"/(f/^Z°,^ ). 
Combining this with the above results, we find that the action (2.2) takes the form 

S^^[Z'- 7,^] = - / dS\l-llh{mn + w) . (2.10) 

JM 

in the Eulerian picture. Here, the mass m and number density n are fixed functions of 
while the interaction energy is a fixed function of hij and Z''. The fleet metric hij is taken 
to be a function of the spacetime metric jais and the variables through Eq. ( |2.7| ). 

Observe that the Eulerian form ( |2.10| ) of the action does not depend on the variable Z^. 
This is a consequence of reparametrization invariance: A reparametrization of the world 
lines induces the transformation 6Z^ = e while leaving the other variables Z* alone. The 
action (2.10) is invariant precisely because it is independent of Z°.[] This gauge freedom is 
removed simply by dropping Z^ from the list of dynamical variables. Thus, we view the 
action (2.10) functional of Z* (and ^ap) only. 



^There is a subtle point here. The action (2.2) is defined for a fixed integration region in S x M; 
that is, for fixed ranges of the integration parameters a and C- The integration region in M is 
determined by the integration region in S x M only if we fix T° at the boundaries. Then, in 
particular, the endpoints in a determine initial and final hypersurfaces in M which we assume to 
be spacelike. Since the range of a in (2.2) is fixed, the action functional (2.10) is defined for the 
class of variables Z^{y) with fixed values on the initial and final hypersurfaces in A4. Therefore the 
gauge freedom in (2.10) consists of variations 5Z^ = e for which e vanishes on the initial and final 
hypersurfaces. In this way we see that the gauge freedom for the actions (2.2) and (2.10) coincide. 
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From the action (2.10), it is straightforward to show that the isentropic perfect fluid 
action S given by Eq. (6.15) of Ref. |jlO|] (or the isentropic case of Eq. (4.20) of Ref. ITll ) 



is equivalent to a "homogeneous and isotropic elastic medium." To do so, consider the case 
in which the mass m is independent of the particle labels Z^, so the particles are identical. 
Also assume that the proper interaction energy density w = w_/ Vh (which is the interaction 
energy per unit proper spatial volume as measured in the rest frame of the matter) depends 
only on the proper particle number density n = n/ \/h (which is the number of particles per 
unit proper spatial volume as measured in the rest frame of the matter). That is, w_ depends 
on and hij only through the combination w_ = \/hw{n/ y/h) for some function 



win) 



The factor {nrn + w)/ ^Jh that appears in the integrand of the action 5*^^ is the proper 
energy density of the medium, which we will denote by p. Our restriction to a homogeneous 
and isotropic medium implies that p only depends on the proper number density, p = p{n) . 
Thus, the action ( p. 101) in the homogeneous, isotropic case may be written as 



SUZ';io.p] = - / Av^p(n) , (2.11) 

where n = n/ y/h. This is precisely the perfect fluid action S given in Refs. [l^jlll for the 
isentropic case, although in those references the number density n was expressed somewhat 
differently. To see that the definition of n in Refs. [[TO| , p!T| agrees with the present definition, 
we use Eq. (2.7) to write the number density explicitly as 



1/2 

(2.12) 



Here, lijk is the antisymmetric tensor density (of weight —1) on S with 6123 = 1. If we 
identify the tensor rjijk{Q of Refs. ||TO| , p!T| with n(C)eijfc, then the expression (|2.12|) exactly 
matches the definition of n given through Eq. (5.1) of [T^ (and through Eqs. (4.23) and (4.26) 
of lll|)- As a result, the ('barred') one component isentropic perfect fluid of Refs. [0,|l3] is 
seen to be a special case of the relativistic elastic medium of DeWitt Q. 



III. CLOCKS AND REFERENCE SYSTEMS 

Reference materials such as elastic media and perfect fluids can be used to provide a 
physical system of coordinates in space. However, by itself, a reference material does not 
provide a complete coordinate system in spacetzme, as all points along a given particle 
worldline are labeled by the same coordinates C*- This can be remedied by adding an 
additional degree of freedom to the particles whose value changes along the worldlines. 
Such a degree of freedom may be called a 'fleet of clocks' and a reference medium coupled 
to a fleet of clocks is said to constitute a reference system. Note that, so far, we have not 
distinguished between 'good clocks' which accurately measure proper time and 'bad clocks' 
whose readings vary along the worldlines in a more complicated way. 



A. Coupling clocks 

The literature contains two different mechanisms for coupling additional 'clock' degrees 
of freedom to a reference material. One of these was used by DeWitt p| and was explicitly 
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described as a coupling of clocks to an elastic medium. The other is implicitly contained 
in the literature on isentropic (single component) perfect fluids |§,10,|ll]|, although the word 



'clock' does not appear in any of these works. Moreover, with a reinterpretation of variables 
and a suitable choice of equation of state, a nonisentropic perfect fluid is actually equivalent 
to a homogeneous, isotropic elastic medium with the DeWitt-type clock coupling. The 
following summary along with the results in Appendix ^ should clarify this situation. 

We first restrict our attention to isentropic perfect fluids. Historically, two different 
action principles were developed for relativistic perfect fluids which both used scalar fields 
( "velocity potentials" ) as the basic variables. These actions were later shown to be equivalent 
pn| . In particular, the action of Ref. [Q, which we will denote by S, was shown to differ 



from the action of Ref. which is denoted by S in Eq. (2.11), by the addition of one degree 
of freedom per space point. This additional degree of freedom is cyclic, so the action S can 



be derived from S by removing the extra degree of freedom through Routh's procedure [|13 

Since the 'un-barred' description of perfect fluids contains an extra degree of freedom per 
space point, and since this degree of freedom changes along the worldlines, the philosophy 
stated at the beginning of this section allows us to interpret this degree of freedom as 
representing a fleet of clocks. For this reason, we refer to the 'un-barred' fluid as a fluid 
coupled to a fleet of clocks. 

For both the clock coupling discussed by DeWitt and the one implicit in the perfect fluid 
literature, the essential idea is to add a pair (9, J) of first order degrees of freedom for each 
particle in the medium. Thus, we consider two fields, 9 (a, () and J(cr, (), and add the first 
order kinetic term 

"da f d^CnJQ (3.1) 



to the Lagrangian picture action ( p.2|) . The Eulerian form of the kinetic term may be 

'/3 



obtained from Eq. (2.9) and the relation 5^ = T"Z*^,^ +T",i Z*,/^ and is given by 



M 



dSnJe,a W^J-^/h . (3.2) 



The clocks are then coupled to the particles by letting either n{() or m{() depend on J in 
the original action. 

To produce the 'un-barred' fluid action of Ref. [^ p!OtpT| , the clocks are coupled by letting 



n(C) Jr{C) in the action (2.11) and adding the kinetic term (3.2). A detailed explanation 
of how this produces the 'un-barred' action of Ref. P JTO|JTl[] is given in Appendix 0. This 



method of clock coupling does not, in general, lead to a 'good' set of clocks. This can be 
seen from the equation of motion for J, which is 

9,„ f/" = p'{Jn/Vh) (3.3) 

where p'{n) = dp(n)/dn. The relation ( |3.3| ) shows that the rate of advance of the clock 
variable 9 relative to the flow of proper time may depend on the particle density n. Thus, 
for clocks coupled in this way, the internal clock mechanism is not shielded from the external 
pressures and forces between particles. For some equations of state p(n), such clocks may 
not even run monotonically in proper time. 

Because the method n —>■ Jn of clock coupling does not in general yield a good set of 
clocks, we will focus on the method described by DeWitt 0. In this method the mass m is 
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allowed to depend on J. Thus, the mass of each particle is no longer a fixed constant and, 
in fact, it acts as a Hamiltonian m{(^, J) that drives the motion of the clock Q{C) attached 
to particle C- the Lagrangian picture, the action Sll[T°', J, 9; 'jais] is obtained by adding 
Eqs. (2.2) and (3.1) and taking the number density, mass, and internal energy density to 
depend on (\ J, and hjk through n = n(C*), ^ = ^{C\ J)^ and w. = w{Ci hjk). The action 
also takes a simple form in the Eulerian picture: 

S,^[Z\ J, 6; 7a;3] = A\/-7//^{li^0,a f/" -{nm + w)] , (3.4) 

where n = n{Z^), m = m{Z\ J), and w = w{Z^, hjk)- Note that the equation of motion for 
J shows that the clock satisfies 



e,„ f/" = dm/dJ , (3.5) 

while the equation of motion for G shows that J (and therefore dm/dJ) is a constant 
along each particle world line. Thus, for the DeWitt-type coupling, increases in direct 
proportion to proper time along the worldline and always defines a good clock. We see that 



( p.3| ) and ( |3.5| ) coincide when the medium is dust, since in that case p' is a constant. 

Also note that many different clock Hamiltonians m{(, J) lead to equivalent results. This 
is because, for any invertible function /(J), the replacement of J by J = f{J) and G by 
G = Q/f'{J) changes the kinetic term (3.1) only by a boundary term but changes the clock 
Hamiltonian from m{(,J) to m{(,f{J)). Thus, any two clock Hamiltonians mi and m2 
related by mi{(, J) = m2(C, f{J)) for invertible / are equivalent. 

Finally, consider a homogeneous and isotropic reference system with DeWitt-type clock 
coupling. As before, the internal energy density has the form w = \^w{n/\/h), and now 
the mass is a function of J only: m = m{J). From Eq. (|3.4|), the action becomes 



S^^[Z\ J, G; 7a;3] = / Av^W0,a U'^ I ^ ' p{n/ J)] , (3.6) 

where pin^ J) = nm{J) + w{n). This action is equivalent to the action (refered to as the 
'hybrid action' in Ref. |TT[]) for a nonisentropic perfect fluid with equation of state p{n, J) = 
nm{J) + w{n). In the perfect fluid literature, the variable J is interpreted as the entropy per 
particle and the variable G is interpreted as the thermasy |12[ . (The thermasy is a variable 
whose gradient along the particle worldlines is proportional to the local temperature.) The 
connection between the action ( ^.6] ) above and the action of Ref. can be established 
easily. One simply uses Eqs. (2.7) and (2.8) to show that the quantity / hnU"' that 
appears in the action (3.6) is the same function of and 7^/3 as the quantity J" that 



appears in the 'hybrid action' of Ref. [11 . 



B. Hamiltonian formulation of reference systems 

We would like to study the diffeomorphism invariance of general relativity by using the 
reference system defined by De Witt's method of coupling clocks to an elastic medium. Since 
we will examine this issue from the canonical perspective, we first perform a Hamiltonian 
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analysis of the reference system. In order to make contact with the usual Hamiltonian 
description of gravity, we will work initially in the Eulerian picture. 

We begin by foliating spacetime A4 with a family of hypersurfaces labeled by the 
parameter t. If x"' are coordinates on a hypersurface, then the spacetime coordinates are 
related to (t, x") through mappings y°' = Y°'(t,x). As usual, the lapse function A^-*-, shift 
vector A^", and spatial metric gab are related to the spacetime metric by 

r° = N^n'' + N'^Y'^.a , (3.7) 

9ab = Y",a7afsY^.l> , (3.8) 

-n'^n^ + Y-^^g-'^YP,, , (3.9) 



a/3 
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where is the unit normal of the hypersurfaces. The quantities A^-*", A^", gab, and n° are 
functions of t and The spacetime metric depends on t and x° through the mapping Y"'; 
that is, ■ja/s = 1af3{Y{t,x)). Note that we are now using the 'dot' to denote d/dt. 

The variables Z'^[y), B, and J are spacetime scalars. They can be pulled back from Ai 
to S X iR by the mapping y"(t,a;) to yield t-dependent scalars on S which we denote by 
the same kernel letter: for example, Z^{y) Z^{t,x). With the definitions above it is not 
difficult to show that (with a slight abuse of notation) 

Z\^n'^ = {Z'-N^Z\a)/N^ , (3.10) 

Z\aY"',a =Z^,a , (3-11) 

Z\'^ = -n''Z\fsn^ + Y'',a9'"'Z\, . (3.12) 

Similar relations hold for B and J. 

The function Z'^{t,x) is a t-dependent mapping from the space S to the matter space 
S. It is useful to consider the inverse mappings X"-{tX), defined by C = Z\t,X(tX))- 
The mappings and X'^ are related by the identities 5* = Z\a^°',j, ^b = ^""nZ^b, 

= — Z*,aX",t, and X°- = —X"',iZ\f The notation used here is somewhat abbreviated 
in that we have omitted any explicit specification of the functional dependences. We will 
continue this practice below. 

Now consider the clock kinetic term (|3.2|) . From Eq. ( p.8|) we have 

e,«f/" = -^e„^^5e,"z^^z^^z^%■fc (3.13) 

which, using Eq. (3.12) and the identity n'^eap^sY'^ ,aY^ ,i,Y^ ,c = eabc, takes the form 

1 .A . 1 



= j^{e'N-Q,a)^e^''Z\aZ\bZ^,e,,k 
1 



^^^,Z' - N-Z^aV'^Q^aZ^ ,bZ\,e,^k . (3.14) 
It is convenient to write this result as 

B,„ = VQ/N^ + Viy - NyN^)e,a , (3.15) 

where 
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r = -n^U'^ = ie'^''^Z^„Z^6Z^,e,,fc , (3.16) 

ya ^ _i_(^a ^ ^a) ^ _ _l_^abc ^^.^^^ Z\c {Z' - Z\a) ■ (3.17) 

Here, V"' is the spatial velocity of the particles as measured by the observers who are at 



rest in and F = l/y 1 — V"-gabV^ is the relativistic 'gamma' factor that characterizes the 
boost between these observers and the co-moving reference frame. 

Using Eqs. (2.7) and (3.16), we find that the fleet metric and its determinant can be 
expressed as 

h'^ = Z\aZ',,{g'''-V''V') , (3.18) 
v^= v^r/(detZ\J . (3.19) 

These results, along with Eq. (3.15), show that the action (3.4) can be written as 

S[Z\J,Q-N^,N\gab] 

= J dt J^d^x{detZ\a)^nJe + nJ{N^V'' - N^)e,a-N^{nm + w)/T^ . (3.20) 

It is now straightforward to proceed with the canonical analysis. The momentum conjugate 
to O is (det Z\a )rJ, which we will denote by 11. The momentum conjugate to is 

dC r 1 

P, = -^ = -{detZ\a) r{nm + w)V, + t,,Z\,Z',,VyT + nJQ,, , (3.21) 

where t^j = 2{dw_/ dW^) is the Lagrangian frame stress tensor 0. Note that appears in 
this expression only in the combination = —{Z^X°',i—N"-)/N-^. Thus, the solution for 
as a function of the canonical variables, lapse, and shift has the form 

Z' = {N" - N^V'')Z\a , (3.22) 

where we view V"' as a function of the canonical variables through equations (3.21) and 
(3.18) for Pi and h'K 

Collecting together the above results, we find the Hamiltonian form of the action for an 
elastic medium coupled to clocks in the Eulerian picture: 

S^^[Z\P,,e,U;N^,N-,gab] = J dt J^d'xi^PiZ' + Ue - N^H"]^ - N^H^^ , (3.23) 

where 

H^ = p,z\,+ne„ , (3.24) 

n'l = (det Z\,)[T{nm + w) + U^Z\a VZ^^ V'' jV] . (3.25) 

Here, F is defined in Eq. (3.16) and the matter spatial velocity V"" and fleet metric /ijj 
are defined as functions of the canonical variables through Eqs. (3.18) and (3.21). The 
clock Hamiltonian, number density, and interaction energy density have dependences m = 
m(H/n(det Z\a ), Z^), n = n{Z'), and w = w{hij, Z^). 
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IV. REFERENCE SYSTEMS IN GENERAL RELATIVITY 



Thus far we have treated the gravitational field as an external background. In this 
section, we shall include gravity as a dynamical field and study our reference system as 
coupled to general relativity. 



A. Coupling clocks to the gravitational field 

Since the action for the matter/clock system contains no derivatives of the spacetime 
metric, the action for that system coupled to gravity is obtained by adding its action to the 
gravitational action. In the Eulerian setting, we have 

= Jdt l^£x!^p,z' + uQ + p'''ga,-N^{n':^ + nl)-N''in': + ni)^ , (4.1) 

where Ti^ and Ti.^ are the familiar constraints of vacuum general relativity |jl4| . The coupled 
system is subject to the constraints 

H^ = n'l + nl = , (4.2) 

na^n:' + ni = o, (4.3) 

since the lapse function iV"*" and shift vector A^*^ are varied in the action principle. 

It is also useful to consider the coupled system in the Lagrangian picture. The Hamil- 
tonian form of the action, namely Shl, can be obtained from a 3+1 decomposition of the 
Lagrangian form of the action in the Lagrangian picture, namely S^i^. Alternatively, Shl can 
be obtained from the Eulerian picture action She by performing a point canonical trans- 
formation from Eulerian to Lagrangian variables. The details of this transformation are 



described in Ref. |Tl|] for the perfect fiuid case and are essentially the same here. Thus, we 
define the new matter variables X'^(C) such that Z^{X{()) = (\ We also define the new 
clock variable 0(C) = ©(-^(0) ^"^^ the new gravitational field variable 



^7.,(C)=x^,(c)x^,(C)^7a6(x(c)). (4.4) 

This is the spatial metric expressed in the Lagrangian system of coordinates induced on the 
slices S by the matter. The new conjugate momenta are defined by ||11||: 



P„(C) = -(detX^,)K(X(C)), (4.5) 

^(c) = (detx^,)^(x(c)) , (4.6) 

fiO = idetX^,, )Z\,{XiO)Z^,,iXiO)p^\X{0) . (4.7) 
With this transformation, the action (4.1) becomes 

SHdx'',Pa,e,Ii,g^,,f,N^,N^] 

= j dtj^ rf=^c{Za^" + ne + fg,, - N^{X)Ul + A^"(^)Za} , (4.8) 
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where the lapse and shift are functions of X°'{() and the Hamiltonian constraint is 2i± = 
(detX°,j )7i_L. In terms of the new canonical variables, we have 

K± = T{nm + w)+ UjZ\a V'Z^^ /V + 2i\ ■ (4-9) 

Here, we have used 21i(C) = (det X",j to denote the pullback [by the mapping 

of the gravitational contribution to the Hamiltonian constraint from E to S. Since 
the Lagrangian gravitational variables gij{C,) and p*-^(C) are obtained from the Eulerian grav- 
itational variables gab{x) and by the pullback mapping from S to S, the term is 
constructed by replacing the Eulerian variables in Ti.^ with the corresponding Lagrangian 
variables. Likewise, we can define the pullback of the gravitational contribution to the 
momentum constraint by Kf (C) = (det T^f (X(C)). 

The velocity V"" and fleet metric that appear in the Hamiltonian constraint (4.9) 
must, of course, be expressed in terms of the new canonical variables. It turns out to be 
convenient to work with which is the velocity expressed in the Lagrangian coordinate 

system, in place of V"". In Appendix ^ we show that hij and VaX°',i can be expressed in 
terms of H, PaX°',i, gij, and p'^^ . 

In the action (4.8), the lapse function and shift vector are functions of X"'{(). We can 
define the new variables A^"'"(C) and N^{C) = ^"•{X^O) by mapping N-^(x) and N"-{x) 
from E to S. The action then reads 

S„,[X'^,Pa,Q,U,g,„f,N^,N'] 

= JdtJ^ d^C^EaX" + ne + p'^gij - N^H^ - N'Hij , (4.10) 

where the Hamiltonian and momentum constraints are 

= T{nm + w) + r¥^KX"« VbX\j +nl , (4.11) 
m = -PaX^n . (4.12) 

Here, the Hamiltonian constraint Ti^ is obtained from Eq. (4.9) and the identity (B2) of 
Appendix 0. The momentum constraint 2ii equals the result obtained by mapping Tia to 
the matter space 5; see Eq. (4.5). 

As a final remark, let us point out that there are two symmetries that play an important 



role in the description of reference systems [p|, p!Ol , pT| . The first is the freedom to reset the 



clocks. That is, the action is invariant under a transformation in which the zero point of 
each clock is changed. The infinitesimal version of this transformation is generated through 
the Poisson brackets by the functional 

Q[^] = J^d^(m = J^d^x^{Z{x))U{x) . (4.13) 

Thus, the clock at C is transformed according to SQ{Q = Qld]} = 'i9(C)- The 

canonical variables in the Lagrangian picture, other than 0(C), remain unchanged under 
this transformation. In the Eulerian picture, we have SQ{x) = 'i?(Z(x)) and SPi{x) = 
—'d,i{Z{x))Il{x), with the other canonical variables unchanged. The second symmetry is 
an invariance with respect to changes in the Lagrangian coordinate labels. That is, the 
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action is invariant under diffeomorphisms of the matter space S. The infinitesimal version 
of this transformation is generated by 

QIC] = Id'CCiEaX'^n+Wn+Kf) = - jj^xC{z{x))Pi{x) , (4.14) 

where the vector field ^*(C) is the infinitesimal generator of the diffeomorphism. The canoni- 
cal variables in the Lagrangian picture just transform according to their 5-tensor character. 
For the Eulerian picture variables we find 5Z\x) = —^\Z{x)) and SPi{x) = {Z{x))Pj{x), 
while the remaining variables (which are tensors on S) are invariant. 

B. New constraints and the gravitational Hamiltonian 

When the Hamiltonian constraint can be resolved for the clock momenta, a new set of 
constraints may be introduced that allows the system to be deparametrized. We begin by 
recalling that, as shown in Appendix B, the fieet metric hij and velocity VaX°',i are functions 
of n, Qij, and 2i.i = ~P^X°'n- It follows that the Hamiltonian constraint (4.11) depends 
on the gravitational canonical variables g^j and p*-^, the clock momentum H, and also the 
particle canonical variables in the combinations given by the momentum constraints Tii- We 
can set 2i.i equal to zero in 2ij^ without changing the content of the constraints ?i± = 0, 
Tij = 0. Then the Hamiltonian constraint Tij^ depends only on gij, and H. 

Let us assume that we can resolve the constraint 2£± = with respect to H. We would 
then obtain a new constraint that has the form 

n^=R+h[g,„p:^] , (4.15) 

where the true (gravitational) Hamiltonian h is a functional of the gravitational variables 
only. The constraints ?i| = and 21i = constitute a complete set of constraints for the 
system which are equivalent to the original Hamiltonian and momentum constraints. We 
can smear 7i| with a prescribed function N\Q on S to form the functional H[N^. The 
smeared constraint H[N^ generates, through the Poisson brackets, changes in the canonical 
variables that result when the hypersurface S is displaced along the particle world lines by 
the clock time N^{Q. 

The constraints 2it h^ive vanishing Poisson brackets among themselves and with the 
momentum constraints 2i.i- As usual, the momentum constraints form a representation of 
the Lie algebra of spatial diffeomorphisms. (Note that, since the gravitational variables gij, 
p^^ and clock variables G, H are independent of S, they have vanishing Poisson brackets 
with 2ii-) Since there are no nonvanishing cross terms in the brackets of with itself, the 
gravitational Hamiltonian h must have vanishing brackets with itself: {h{() , h{(')} = 0. We 
should emphasize that this property holds for any gravitational Hamiltonian h that can be 
derived through the coupling to an elastic medium with clocks; that is, for any choice of 
internal energy density w = w{(\ hjk). In practice, the equations that determine h are not 
solvable by analytical means except for the simplest of cases, such as dust {w = 0).| 



^Some other matter couplings also generate true Hamiltonians that depend only on the gravita- 
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Finally, note that the constraint (4.15) allows the system to be deparametrized. That 
is, consider the action S^l with N^T^x replaced by A^^Hf- The particle variables X° and 

are completely decoupled from the gravitational and clock variables and can be dropped 
from the action. The variables N'^ and n can be eliminated through the solution of their 
equations of motion, and the parameter t can be replaced by the clock time 9 = 0(t). This 
leads to the action 

S[9^,,f] =JdeJ^ d'c{fg., - h[g,,,f]} , (4.16) 

where the dot denotes d/dO. This action for the gravitational field contains no gauge (dif- 
feomorphism) invariance. 



C. The case of dust 

We will now explicitly display the results for dust, where the interaction energy density 
vanishes (i£ = 0). In this case, as in Ref. |p, the Hamiltonian constraint (4.11) becomes 



rt^ = ^{nmf + VSg'^Wj + K\ ■ (4. 17) 

Recall that n = n{() and m = m{(,Il{Q/n{Q). Let us consider various possible choices 
for nm. With the clock Hamiltonian m(C, J) = J/k, where /c is a positive constant, then 
nm = n/ k and the new constraint takes the form 2i| = H + h with the gravitational 
Hamiltonian 



h = Tk^{niy-nfgmh (4.18) 

As described in Ref. , the presence of the square root creates a serious difficulty for defining 
a quantum theory. In particular, it implies that the physical Hilbert space must be restricted 
to those states for which the operator (^x)^ ~ itiQ^^itj is nonnegative, and the physical 
observables must be restricted to those operators that keep the states in the physical Hilbert 
space. The difficulty is that the obvious candidates for physical observables, the metric 
9ij ~ 9ij^ conjugate p*-' = —i6/6gij acting by multiplication and differentiation, do 

not satisfy this criterion. 

While different choices for the clock Hamiltonian are classically equivalent, as long as they 
are defined for the same range of J, they do not necessarily lead to the same quantum theory. 



Thus, we now investigate other possibilities. With the clock Hamiltonian m{(, J) = J J/k 



(where we assume J > 0), we find that nm = J nli /k and the gravitational Hamiltonian is 



h 



mf-mg'm ir. (4.19) 



This result suffers from the same difficulty as the square root Hamiltonian (4.18) if we 
insist that the momentum H should be interpreted in terms of a real, nonnegative clock 



tional variables. Examples that can be solved analytically include the massless scalar field and 
certain perfect fluids with 'bad' clock coupling ||l6|| . 
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Hamiltonian m. In that case n must be nonnegative (assuming the number density n is 
positive), and the new constraint 2i.i = H + h = imphes that h of Eq. (4.19) must be 
nonpositive; we again have the requirement (21_l) — TlLig^''TlLj > on the quantum theory. 
On the other hand, it is not obvious that the condition n > must be kept, so one might 
argue that it should be dropped. This leads to a quantum theory that is free from the 
difficulties of the condition (Kj,)^ ~ 'tLiQ^^'tLj > 0, but there might arise a problem of 
interpretation for the theory since the variables and H would no longer have a simple 
interpretation in terms of a real clock Hamiltonian. 

Clearly, by choosing the clock Hamiltonian m{(, J) appropriately, we can arrange for 
the gravitational Hamiltonian h to be given by any invertible function of (211)^ ~ 2M9^'''t^j 
(with n appearing as necessary to keep the density weights balanced). In particular, for the 
clock Hamiltonian m((^, J) = (J/kY^^ (again with J > 0), we have nm = rv^^^ill/ky^^ and 



h = -k 



2 



mr-ma'm /{nf. (4.20) 



This choice for the gravitational Hamiltonian appears to avoid the problems encountered 
in the other cases. In particular, the clock momentum remains nonnegative, H > 0, for 
all solutions of the constraint Ti.] = 0, even on states for which {H±) ~ tLig^-'tLj is a 
negative operator. However, some difficulties with interpretation do arise if we attempt to 
take seriously the constraint in its original form 2i± = 0. In that case, for states satisfying 



the constraints and on which (2i±)^ — iLiQ^^lLj is negative, we find that rn? < 0. Thus, for 
such states the clocks are tachyonic.0 One might also be concerned that there is something 
pathological about a Hamiltonian h that is eighth order in the gravitational momenta. 
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APPENDIX A: PERFECT FLUIDS AS FLUIDS COUPLED TO CLOCKS 

By making the substitution n Jn in the action (2.11) and adding the kinetic term 
(3.2), we obtain an action for an isentropic perfect fluid coupled to clocks: 

S[Z\J,Q-^^^]= I dS\-V^p{Jn/^)+nJQ,^U''J^^\. (Al) 

Here h and U"" are expressed in terms of the label flelds through Eqs. ( |2.7| ) and (^.8p, 
respectively. In this appendix we show that this action is equivalent to the (isentropic) 'un- 
barred' perfect fluid action of Ref . P JT(]||TT[ ] . In notation consistent with the present work, 
that action is 



^Said differently, if the mass is still defined by the real branch of (J/fc)^/^, it will not satisfy 



2i_|_ = for all solutions of the new constraint 2i| = 0. 
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M 



d'y -v^p(|ir|/v^) + ir"(e,,+AiZ^ 



(A2) 



where \K\ = \/—K°'Ka is the norm of the timehke, future directed, spacetime four-vector 
density i^". (Here, our i^", Aj, and correspond respectively to J'*, and of Ref. 
and to J", -Wk, and Z*^ of Ref. 

The essential step in comparing the actions (Al) and (A2) is to note that the equations of 
motion for some of the fields in the action (A. 2) can be solved algebraically for those fields 
as unique functions of the other fields. These 'dependent' fields contain only redundant 
dynamical information and, when the corresponding equations of motion are solved and the 
solutions inserted into the action (|A2|) , the resulting action is equivalent to (|A2|) . We will 
see that this procedure produces the action (|A1|). 

In order to identify the proper equations of motion to solve we must first change coor- 
dinates on field space by writing i^" = /tf/", where f/" is a unit future pointing timelike 
four-vector. Thus, the action (A2) can be written as 



M 



-7p(«:/v^)+«:f/°(e,„+A,Z\„ 



(A3) 



The variations in If^ must preserve the normalization condition W^Ua = — 1. 

Note in particular the equations of motion that follow from varying Aj and ?7": 

6S 
SS 








f/"Z\« = 



(A4) 

(A5) 
These 



Here, C is a proportionality constant that arises due to the restriction on SW^ 

equations are easily solved for Aj and The solution for is just the expression 
The solution for Aj is found by contracting Eq. (A5) with and using Eq. (A4). The 

result is 



A,; 



(A6) 



where the fleet metric is defined in terms of by Eq. (2.7). The solution for [/" and Aj can 
be substituted into the action (A3) to yield the equivalent action 



M 



d'y 



(AT) 



where is the functional of given by Eq. (p.8|) . If we now identify k = Jny—'y/h 
(another change of coordinates on field space), this is just the action (|Al|) . We see that the 
'un-barred' fluid of Ref. P, p!0| , [rT| may be interpreted as an isentropic perfect fluid coupled 
to a fleet of clocks. 



APPENDIX B: FLEET METRIC AND MATTER VELOCITY AS FUNCTIONS 

OF THE CANONICAL VARIABLES 

The velocity VaX°-,i and fleet metric hij that appear in the Hamiltonian constraint (4.11) 
are expressed in terms of the Lagrangian picture canonical variables as follows. First, com- 
bine Eq. (3.21) for Pj with Eq. (4.5) for to obtain 
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Pa = r(rim + w)Va + t,jZ\a ,b V'/r - (det X\j )ni . (Bl) 

By solving Eq. (3.18) for g"-'^ and inserting the result into Z\a,V°' = Z\ag°''^Vh, we obtain 
the identity 

Z^aV ^r^h'^VaX^j . (B2) 

Then Eq. (Bl) becomes 

Z„^"n + Kf = r [{nm + w)Sl + t,,h''] VaX\j , (B3) 

where TYf is the puUback to S of the gravitational contribution to the momentum constraint. 
With the identity (B2), one can confirm that 

% = + ^2(KX^, )(H^',i ) (B4) 

is indeed the fleet metric; that is, hij is the inverse of the inverse fleet metric (3.18). 

Now observe that the gamma factor F = 1/ ^/l — can be expressed in terms of VaX"',i 
and h'^K This can be seen by again solving Eq. (3.10) for g"-'^ and inserting the result into 
y2 = Vag^'^Vb. This yields V^/T^ = {V''X'',i)h'^ {VbX\j), which can be solved for T as a 
function of VaX°-,i and k^^ . We therefore see that, in principle, Eqs. (B3) and (B4) can 
be solved for h^j and V"'X"-,i as functions of the canonical variables. In particular, hij and 
V^X",j depend on n [which is contained in the argument of the clock Hamiltonian m in 
Eq. (B3)], P_aX°',i [which appears on the left-hand side of Eq. (B3)], gij [which appears 
explicitly in Eq. (B4) and on the left-hand side of Eq. (B3) in the combination 2if ] and p^^ 
[which appears on the left-hand side of Eq. (B3) in the combination Hf]. 
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